First-principles study of phenyl ethylene oligomers as current-switch 
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We use a self- consistent method to study the distinct current- switch of 2 -amino-4-ethynylphenyl- 
4 / -ethynylphenyl-5 / -nitro-l-benzenethiol, from the first-principles calculations. The numerical re- 
sults are in accord with the early experiment [Reed et al., Sci. Am. 282, 86 (2000)]. To further 
investigate the transport mechanism, we calculate the switching behavior of p-terphenyl with the 
rotations of the middle ring as well. We also study the effect of hydrogen atom substituting one 
ending sulfur atom on the transport and find that the asymmetry of TV curves appears and the 
switch effect still lies in both the positive and negative bias range. 

PACS numbers: 73.23.-b, 85.65.+h, 31.15.Ar 



As Moore's law predicted, the density of transistors 
on integrated circuits doubles approximately every 18 
months, and the traditional semiconductor devices will 
approach to their limit. In this century, the develop- 
ment of electronics may enter a new era of molecular 
devices. Compared with the traditional solid electronic 
devices, the molecular devices have many advantages. 
The molecular chips will be smaller than the silicon chips 
in three orders and their operation speed will increase 
evidently without increases in price. The concept of 
the organic molecular electronic devices was first pro- 
posed over a quarter century ago by Aviram and Rat- 
ner ^J. However, the molecular devices did not attract 
much attention until the measurement of the electronic 
conduction through the single phenyl dithiol (PDT) by 
Reed et al. in 1997 |2|. From then on, the search for 
new active molecular devices becomes a worldwide effort 
H 0, 0, 0, El El E 13 El E3 Molecules are 
intrinsically different from traditional devices, and one 
expects that some of their transport properties would 
be unique. An attractive feature of organic molecules 
for application to electronics is the possibility to con- 
trol electronic transport by using electrical or photic sig- 
nal. Several experimental groups reported some molec- 
ular switch devices 0, 0, ^| , which are mainly in- 
cluded into two categories. One is the electronic switch 
devices which control the conversion between state "0" 
and state "1" by the reversible field or the current pulse. 
The other is the photic switch devices which achieve the 
open or close states of current by the laser field. Reed et 
al. E3 used a molecular switch consisting of three aro- 
matic phenyl rings in series. The two hydrogen atoms 
of the middle ring are substituted by acceptor group 
NO2 and donor group NH2, while the whole molecule is 
chemisorbed onto the contact surfaces of gold leads. By 
adjusting the gate voltage, the current can be controlled 
due to rotation of the middle ring with respect to the 
side rings. The corresponding theoretical works investi- 
gated the similar switch system [H 0, However, 
in these works, the macroscopical electrodes were not 



treated on an equal footing with the molecular devices, 
and the whole open system, lead-molecule-lead, were not 
treated self-consistent ly with the charge effect included. 
In the quantum system, the charge effect is crucial to 
the transport property. The rigorous treatment of the 
molecular device calls for combining the theory of quan- 
tum transport with the first-principles calculations of the 
electronic structure. In this article, we use density func- 
tional theory (DFT) and nonequilibrium Green's func- 
tion to study the current-switch behavior of the phenyl 
ethylene oligomers. 

For the lead-molecule-lead system, coupling between 
the molecule and the leads plays a crucial role in quantum 
transport. Being computationally tractable, the whole 
system is partitioned into the molecule part and the lead 
part, so that these two parts are dealt with separately. 
The nonequilibrium Green's function theory provides a 
powerful method to give a full description of transport 
phenomena. The molecular Green's function is expressed 
as follows 
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where G Ml Sm, Fm are the retarded Green's function, 
overlap matrix and Fock matrix of the molecule part re- 
spectively, the retarded self-energy of the left 
(right) electrode, is calculated from the surface Green's 
function (SGF) g* (g*) 
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with i = 1,2. The coupling matrices Sm% and Fm,% 
are extracted from the DFT calculation for the extended 
molecule (molecule with 3 Au atoms on each side). 

The density matrix of the open system is the essential 
function of the whole self-consistent scheme. It can be 
achieved by the Keldysh Green's function 
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with the advanced Green's function G A = {G R )\ the 
broadening function of the left (right) lead Ti (r 2 ). 
The Fermi distribution function of the left (right) lead 
h (/ 2 ) is expressed U(E) = i/( e (E-^)/kT + ^ with 
111 = Ef — ^eV, /12 — Ef + \eV we regulate Z axis 
pointing from the left to the right, so V > (V < 0) de- 
notes the longitudinal electrical field direction along the 
positive direction of Z axis (along the negative direction 
of Z axis). But, this expression is only right for the case 
that the molecule is symmetrical to the two contacts. If 
not symmetrical, the voltage drop is also not symmetri- 
cal to the two interfaces. Fermi level of the bulk Au is 
Ef. In our work, Ef is -5.1 eV which is adjusted around 
its work function (5.31 eV). Since the relative dielectric 
constant of the molecule is much larger than 1, the volt- 
age drops across the interface between gold atoms and 
the sulfur atom, while the electronic potential is almost 
flat. So, the potential zero point can be set at the center 
of the molecular device. 

After obtaining the density matrix self-consistently, we 
calculate the transmission function T(E, V) of coherent 
transport 

T(E, V) = Tr(r 1 G i? r 2 G A ). (5) 

In order to achieve this goal, we extended the inner loop 
in the standard quantum chemistry software Gaussian03 
[2^ to the loop containing the lead-molecule-lead open 
system under bias. The self-consistent procedure starts 
from a guess for the density matrix of the open system, 
which may be obtained from the converged density ma- 
trix of Gaussian03 calculation for the isolated molecule. 
We feedback the density matrix to the Gaussian's main 
program as a subroutine to obtain the new density ma- 
trix. The iterations continue until the density matrix 
converges to the acceptable accuracy. Then the den- 
sity matrix is used to evaluate the terminal current and 
density of states (DOS) of the open system [2J,[2^. In 
the calculations, we adopt DFT with B3PW91 exchange- 
correlation potential and LANL2DZ basis to evaluate the 
electronic structure and the Fock matrix. The basis set 
associates with the effective core potential (ECP), which 
are specially suited for the fifth-row (Cs-Au) elements 
with including the Darwin relativistic effect. 

In Reed's experiment 0, the middle ring of the 
molecule 2 -amino-4-ethynylphenyl-4 / -ethynylphenyl-5 / - 
nitro- 1-benzenethiol can rotate with respect to the side 
benzene rings under control of the gate voltage. The 
dipole composed of the acceptor NO2 and the donor NH2, 
which are attached to the middle benzene, is driven by 
the external field. Our DFT calculation at LANL2DZ 
level achieves the optimized molecular structure, where 
three aromatic phenyl rings are almost in the same plane, 
the carbon bond length of the alkyne in the middle of two 
rings is 1.24 A, the angle of O-N-0 is 128° and the an- 
gle of H-N-H is 118°. For the sulfur atoms chemisorbed 
onto the gold leads, there are two kinds of the geometric 



configurations. The sulfur atom sits directly on the top 
position of the surface gold atom or the hollow position of 
three nearest-neighbor surface gold atoms. As a conven- 
tionally accepted picture the latter is adopted in this arti- 
cle. The perpendicular distance between the sulfur atom 
and the Au FCC (111) surface plane is 2.0 A, an usually 
acceptable distance. Due to the precalculation, the tem- 
perature effect is not distinct for the short molecule, so, 
we assume zero temperature in our calculation for sim- 
plicity. It is emphasized that since our system is a closed 
shell system (the number of a electrons is equal to that 
of /3 electrons), we adopt restricted calculation which en- 
sures that p a = in order to spare calculation time. 
We investigate the molecules under three configurations 
corresponding to the rotation angle a = 0°, 30°, and 60° 
respectively. DOS and T corresponding to the different 
values of a are given in Fig. 1, which accounts for the 
mechanism of the molecular switch. Both HOMO and 
LUMO are shifted and broadened with the molecule at- 
tached by the metallic leads, so that Ef is almost in the 
middle of HOMO-LUMO gap, a little closer to HOMO, 
which means electrons will be responsible for the initial 
rise of the current. The electron incident from the left 
lead with the energy meeting the one of molecular states 
will enter the resonant channel. The states with small 
DOS in the HOMO-LUMO gap is called as the metal- 
induced gap states (MIGS) arising from the hybridiza- 
tion of gold surface states and molecular HOMO, LUMO 
states. They offer some electrons (about 0.32 electrons) 
to the contacts, which causes the molecule to be a system 
with positive charges. When bias is applied, the system 
in the nonequilibrium state adjusts its energy levels by 
the charge effect to prevent the loss of electrons. 

The transport characteristics of the molecule-lead sys- 
tem with the rotated configurations are obtained from 
the DOS curves. The singular points of the Green's func- 
tion, which are obtained from equation + + £2)6* = 
SCX, show the position of the levels of the open system. 
In the case of a = 0°, HOMO is -6.53 eV, LUMO is - 
3.60 eV, and the HOMO-LUMO gap is 2.93 eV. In the 
case of 30°, HOMO is -6.54 eV, LUMO is -3.56 eV, and 
the HOMO-LUMO gap is 2.98 eV. And for 60°, HOMO, 
LUMO and the HOMO-LUMO gap are -6.73 eV, -3.51 
eV, and 3.22 eV, respectively. With the increase in the ro- 
tation angle, HOMO and LUMO of the open system keep 
away from Ef step by step, and the transmission func- 
tion value at the Fermi level gets smaller. The transmis- 
sion functions are obtained for the three configurations, 
T(E f ) = 0.19 xlO- 2 , 0.14 xlO- 2 , 0.34 xl0~ 3 correspond- 
ing to 0°, 30° and 60°. The big rotation angle (e.g. 60°) 
results in the heavy decrease in the current (the 90° rota- 
tion causes the switch into the off state, not shown in the 
figure) which demonstrates the switchable function of the 
molecular device. With the molecular device attached by 
the gold leads, the energy levels are broadened and de- 
scend to keep the system electrical-neutral. The electron 
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FIG. 2: I-V curves of 2 -amino-4-ethynylphenyl-4 / - 
ethynylphenyl-S'-nitro-l-benzenethiol with gold contacts 
through sulfur atoms from both sides corresponding to 
a = 0°, 30° and 60°. 
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FIG. 1: 2 -amino-4-ethynylphenyl-4 / -ethynylphenyl-5 / -nitro- 
1-benzenethiol is chemisorbed onto the gold leads through 
sulfur atoms from both sides. DOS (solid) and T (dashed) as 
functions of energy in equilibrium. The vertical line denotes 
the position of Fermi level, (a) rotation angle a = 0°, (b) 
a = 30°, (c) a = 60°. 



number inside the device decreases from 186 (in the iso- 
lated molecule) to 185.6. The basic reason is illustrated 
by the inset of Fig. 2. The inset shows the energy levels 
of the isolated molecule corresponding to the three cases. 
LUMO of the isolated molecule is below the Fermi level, 
and LUMO+1 is above the Fermi level. With the increase 
of the rotation angle, LUMO-LUMO+1 gap gets bigger, 
leading to the same tendency shown in Fig. 1. With leads 
attached, the charge effect makes the energy levels of the 
isolated molecule descend and broadened. The LUMO- 
LUMO+1 gap evolves into the HOMO-LUMO gap of the 
open system and determines the transport capacity. In 
Fig. 2 the I-V curves are steep for a < 30°, since both 
HOMO and LUMO are responsible for the rising of cur- 
rent. The negative differential resistance (NDR) appears 
with the added electrons, when bias gets across -3.6 V. 
From 30° to 60°, the switch effect due to the rotation is 
more apparent than the case a < 30°. 



The p-terphenyl is a good sample to investigate the 
molecular switch properties. Fig. 3 gives the I-V curves 
of p-terphenyl with the anti-clockwise rotations of the 
middle ring with respect to the side benzenes. In the fig- 
ure there are five current curves corresponding to rota- 
tion angle a = 0°, 30°, 45°, 60°, and 90°. The molecular 
device presents an active switch function, although it is 
not easy to make it work in practice. The total energies 
of the isolated molecule corresponding to these angles 
are -712.967, -712.962, -712.952, -712.938, and -712.924 
Hartree, respectively. The energy increases with an in- 
crease in the rotation angle, while the molecule with the 
zero angle is the most stable one. The figure illustrates 
a distinct switch function. From 0° to 90°, the overlap 
between ir electron clouds decreases with the increase 
in the rotation angle, which leads the current reduction. 
For the 0° case, the current of p-terphenyl under the bias 
3.0 V is smaller than the one of single PDT, while it is 
bigger than the latter one above the bias 3.0 V. The ab- 
normal phenomenon comes from the fact that HOMO of 
p-terphenyl is closer to Ef. The current has a big rise at 
the point where the voltage gets across HOMO level. The 
inset shows the molecular orbital pictures of HOMO for 
rotation angle 0° and 90°. Here, the extended molecule 
(3 Au atoms on each side are included) is adopted to 
simulate p-terphenyl connected to gold contacts and its 
HOMO level has the main contribution to transport in 
our bias range, while LUMO does not participate in. It 
is obvious that the HOMO is a delocalized state for 0° 
and it becomes a localized one when the rotation angle 
approaches to 90°. 

Typically the aromatic phenyl molecule is chemisorbed 
on the gold surface through a sulfur atom. J. Chen et 
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FIG. 3: I-V curves of extended p-terphenyl compound cor- 
respond to the different rotation angles. The negative bias 
means the electrical field is along the direction opposite to 
the Z axis. 



al. |26| once investigated the system with a hydrogen 
atom to replace the sulfur atom, and obtained a differ- 
ent transport characteristics, since the electronegativity 
of hydrogen atom is smaller than that of the sulfur atom, 
which makes the coupling between a hydrogen atom and 
one gold contact weaker than that between a sulfur atom 
and the other gold contact. The case is attractive be- 
cause of the asymmetry of I-V curves [27|. We add a 
hydrogen atom on the sulfur atom at the right end to 
make the whole molecule a closed shell system to spare 
the calculation time, without the intrinsical influence on 
I-V curve. The distance between hydrogen atom at the 
left end and gold contact is a little longer than that be- 
tween sulfur atom and gold contact. But this distance 
is not clear. We fix this distance at 2.15 A. The dis- 
tance between sulfur atom and the other gold contact 
is still fixed at 2.0 A. So, the electrostatic potential dif- 
ference V is divided between the two junctions from the 
coupling distances: 52% across the H-Au junction and 
48% across the S-Au junction. The inset of Fig. 4 gives 
the energy levels of the isolated molecule with the re- 
lated HOMO-LUMO gaps: 3.15 eV, 3.23 eV and 3.43 
eV for a = 0°, 30° and 60°. The HOMO-LUMO gap 
increases with the increase in rotation angle, which dom- 
inates the transport behavior. Fig. 5 describes DOS and 
T in equilibrium with the different rotation angle. With 
the gold leads attached to the molecular device, there 
are about 181.2 electrons in the device area. Different 
from the pure S-Au contacts, only LUMO is responsi- 
ble for the initial rise of the current. With contacts, for 
the case of 0°, HOMO is -8.03 eV, LUMO is -4.73 eV, 
and HOMO-LUMO gap is 3.30 eV; for 30°, HOMO is 
-8.06 eV, LUMO is -4.68 eV, and HOMO-LUMO gap is 



3.38 eV; for 60°, HOMO, LUMO and HOMO-LUMO gap 
are -8.08 eV, -4.62 eV, and 3.46 eV, respectively. The 
corresponding I-V curves are shown in Fig. 4 for bias 
region -1.4 V< V < 3.4 V. The distinct switch effect ex- 
ists in the small negative bias, while the electron number 
descends with the increase in bias. The similar switch 
property occurs in the positive high bias (above 2.0 V) 
with the electron number variation from 181.2 to 181.6. 
Because the molecule is not symmetrically connected to 
the gold contacts, the negative bias is easier to get across 
LUMO than the positive one and for the negative bias 
the charge effect does not occur since the left chemistry 
potential is higher than the right one. Then, LUMO ab- 
sorbs electrons from the left side, and emits electrons to 
the right side. The left coupling is weaker than the right 
one, which keeps LUMO normally unoccupied and re- 
sults in the descendence of electron number, without the 
charge effect. So, the current has the initial rise for the 
small negative bias. For the positive bias, the left chem- 
istry potential is lower than that of the right one, which 
makes LUMO absorb electrons from the right side and 
emit electrons to the left one. The abnormally occupied 
LUMO due to the nonsymmetric coupling causes the rise 
of electron number. So, the LUMO is kept unoccupied by 
the charge effect, which results in very small current in 
the small positive bias range (V < 2.0 V). With high bias 
added, the resonant transmission occurred. This process 
is shown in Fig. 6. With the increase of positive bias 
from 1.0 V to 3.4 V, the peaks of DOS and T shift to 
the high energy side continuously. At 3.4 V, the current 
corresponding to 0° has a sudden rise since the peaks' 
shift makes both HOMO and LUMO responsible for the 
transport at this condition. The sudden rise makes the 
switch effect more apparent in comparison with the case 
of 60° rotation, where only LUMO is responsible for the 
rise of the current. Since the distance between the hy- 
drogen atom and the gold contact is not clear, we also 
calculate the DOS and T of the molecule in equilibrium 
with the H-Au distance 2.3 A. Corresponding to 0°, 30°, 
60°, the HOMO-LUMO gaps are 3.31 eV, 3.34 eV, and 
3.43 eV respectively. The small change of H-Au distance 
has little effects on the positions of peaks of DOS and T, 
only causes current smaller. In fact, even this distance 
is fixed at 2.9 A, the switch effect and the asymmetry 
of I-V curves still exist. Our calculation shows that the 
asymmetry effect and switch effect are both stable corre- 
sponding to the different H-Au distances. We can draw 
the conclusion: for the LUMO-based molecule, the cur- 
rent is lower for positive bias on the weak contact, while 
for the HOMO-based molecule, the current is lower for 
positive bias on the strong contact. The conclusions from 
the first-principles calculation accord with the previous 
results obtained by the extended Huckel theory (EHT) 
[2^ . This asymmetry in I-V curves may be used to pro- 
vide a rectifier. 

In summary, we use the self-consistent method based 
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FIG. 4: I-V curves of 2 -amino-4-ethynylphenyl-4 - 
ethynylphenyl-S'-nitro-l-benzenethiol with gold contacts 
corresponding to a = 0°, 30° and 60°. The atoms that 
chemisorbed to the gold surfaces are hydrogen atom and 
sulfur atom respectively. The distance between the ending 
hydrogen atom and gold contact is 2.15 A. 



on the DFT and the non-equilibrium Green's function 
to simulate molecular transport. We solve the trans- 
port problem from the first-principles theory in order 
to predict the transport characteristics of some molecu- 
lar devices and to identify the experimental results from 
the theoretical point of view. In terms of Gaussian03, 
the electronic structures for the molecular device and 
the macroscopic leads are calculated on an equal foot- 
ing, at the same time, the self-consistent iteration cycle 
is extended from the local molecule to the open system, 
by inserting the density matrix of the open system as 
a subroutine. The self-consistent iteration is run until 
the density matrix converges within an acceptable accu- 
racy. We investigate the molecular switch reported in 
the experiment 0, and prove the distinct switch func- 
tion, theoretically. With one sulfur atom replaced by 
the hydrogen atom, the switch effect exists in both the 
low negative bias range and the high positive bias range. 
The calculation result shows that the asymmetry of I-V 
curves can be used to make the rectifier. Further study 
of the influence of gate voltage on molecular transport 
characteristics is in process. 
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FIG. 5: DOS (solid) and T (dashed) as functions of en- 
ergy of 2 -amino-4-ethynylphenyl-4 / -ethynylphenyl-5 / -nitro- 
1-benzenethiol with gold contacts in equilibrium. The atoms 
that chemisorbed to the gold surfaces are hydrogen atom and 
sulfur atom respectively. The distance between the ending 
hydrogen atom and gold contact is 2.15 A. The vertical line 
denotes the position of Fermi level. The rotation angle a — 0° 
(a), 30° (b), and 60° (c). 



facilities. 
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